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A local and renormalizable version of a modified PQCD introduced in previous works is presented. 
The construction indicates that it could be equivalent to massless QCD. The case in which only 
quark condensate effects are retained is discussed in more detail. Then, the appearing auxiliary 
fermion fields can be integrated leading to a theory with the action of massless QCD, to which one 
C^l ' local and gauge invariant Lagrangian term for each quark flavour is added. These terms are defined 

by two gluon and two quark fields, in a form curiously not harming power counting renormalizability. 
The gluon self-energy is evaluated in second order in the gauge coupling and all orders in the new 
quark couplings, and the result became transversal as required by the gauge invariance. The vacuum 
energy was calculated in the two loop approximation and also became gauge parameter independent. 
The possibilities that higher loop contributions to the vacuum energy allow the generation of a quark 
mass hierarchy as a flavour symmetry breaking effect are discussed. However, the decision on this 
' issue needs the evaluation of more than two loop contributions, in which more than one type of 

quark loops start appearing, possibly leading to interference effects in the vacuum energy. 

^ • PACS numbers: 12.38.Aw;12.38.Bx;12.38.Cy;14.65.Ha 

D I. INTRODUCTION 

The understanding of the hierarchy of quark masses is a fundamental open problem of High Energy Physics. In 
t-H ■ the investigation of this question, the theory of particle condensation in field theory had been a basic framework 
p]4l3|. In an effort to explore this issue, an alternative to the standard PQCD, including the presence of quark and 
gluon condensates in the free vacuum state generating the Wick expansion, has been considered in Refs. (l5rl25| . 
The exploration was in great extent motivated by the consideration about that massless QCD convey the strongest 
forces in Nature, as well as the free theory is a massless highly degenerate one for both quarks and gluons. This fact 
rises question about what could be the real intensity of the dimensional transmutation effect [ll|. The first step 
was to search for simple states of the free theory in which large numbers of zero momentum quark and gluo n states 
were created in the standard free vacuum state, on which, then, adiabatically connect the interaction [l6l [l7l l2lj |. 
In the start, the objective was to determine a modification of the standard Feynman rules of QCD embodying new 
• • ■ condensate effects, with the expectation of to attain a theoretical prediction of superconductivity kind of properties 
. ^ ' of the particle mass spectrum underlined in Refs. 0, B|. 

In Refs. [HHpl], some indications about the possible dynamic generation of quark and gluon condensates were 
obtained. In particular in Ref. 22] it was restricted the study to consider only the existence of quark condensates, 
5^ ! with the aim of exploring the generation of large values for them. The motivation was the suspicion about that 
this mechanism might be able in generating a variant of the top quark model, as an effective action for massless 
QCD. In this case, the generating functional Z of the system introduced in Ref. [20] was transformed to an alternative 
representation in which all the effects of the condensates were incorporated in a new vertex showing two quark and two 
gluon legs. This representation implemented the dimensional transmutation effect produced by the modification done 
in the free vacuum state. The results suggested a technical path that could evidence a possible strong instability of 
massless QCD under the generation of fermion condensates. This outcome could furnish an explanation of the particle 
mass hierarchy in a kind of generalized Nambu-Jona Lasinio dynamical symmetry mechanism [H 13, [l2| . However, 
a drawback of the diagram expansion introduced in Ref. [221 ] was that the new vertex, although being covariant was 
a nonlocal one. That is, it was not associated to the four appearing fields defined at the same space-time point. 
Moreover, the evaluated two loop corrections to the vacuum energy turned out to be unbounded from below as a 
function of the quark condensate parameter, a fact that shifted the answer of the question about the prediction of 
dynamical mass generation to higher loop evaluations. The undesirable properties of the new vertex, can be traced 
back to the form of initial state of the free theory employed to connect the gauge interaction. This wavefunction was 
a sort squeezed state constructed with nearly zero momentum quark creations operators acting on the standard free 
vacuum. It has a structure similar to the BCS state in the usual superconductivity theory. The vertex was derived 
from a particular representation in which the Feynman gauge was explicitly employed and the quark free particle states 
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used in the construction were particular zero momentum states in this special gauge. Therefore, the gauge invariance 
of the description was directly broken by the derivation. Henceforth, it came to the mind the possibility that the 
connection of the interaction on a specially defined in a more sophisticated way vacuum state could implement the 
locality and gauge invariance of the functional integral action in the ending theory. In particular coordinate dependent 
states in this same gauge exist that can be imagined to allow the mentioned construction. 

In this work we start form a the results for the generating functional for the massless QCD in which quark and 
gluon condensate parameters were introduced in the free vacuum, for afterwards construct the Wick expansion in Ref. 
[20(. The resulting scheme included integrals over auxiliary boson and fermion parameters which appeared in the 
process of linearly representing the quadratic forms in the sources introduced by the inclusion of the gluon and quark 
condensates. Then, it is firstly observed that the simple promotion of the gluon and quark condensate parameters 
to be space dependent functions makes the action defining the generating functional a local and gauge invariant 
one. Further, arguments are presented suggesting that the connection of the interaction in some sort of more general 
and yet undetermined initial modified vacuum state, could show the functional to be equivalent to massless QCD, 
after the adiabatic connection of the interaction. Afterwards, the discussion in the work is restricted to the case of 
retaining only the fermion condensate auxiliary functions. The study of the case of the gluon condensate functions 
has also appreciable interest in itself for to explore its possibilities in the description of confinement and low energy 
Physics. However, the issue of the quark mass generation possibility is closer to the fermion condensate parameters 
and we decided to consider it first here. In this case, the Gaussian dependence of the quark auxiliary functions can 
be integrated out to give a theory action given by the massless QCD one plus six terms (one for each quark flavour). 
They implement the above mentioned local and gauge completions of the non local and gauge non invariant action 
obtained in Ref. [22| . 

The strongest support we had found for the physical relevance of these new terms simply comes form their structure; 
it seems that they are allowed counterterms of the massless QCD Lagrangian. In other words, if not obstructed by a 
subtle point not noticed by us, it seems that such a contributions can appears in the effective action of massless QCD 
after finite renormalizations of the four legs 1PI vertex determining the two gluon and two fermion fields terms in 
this expansion. Curiously, the new terms does not obstacle power counting renormalizability. This apparently strange 
property can be understood by taking into account that the new Lagrangian term also includes a contribution to the 
quark free propagators, which acquire a more converging behaviour 1/p 2 at large momenta. Thus, the theory remains 
being power counting renormalizable. The above mentioned properties furnish an intrinsic interest to the theory 
resulting after simply adding the new local and renormalizable terms to the massless QCD action, as an alternative 
to the standard QCD. 

Further, the Feynman diagrams associated to the new vertices and quark propagators are defined. As before noticed 
masses for the free fermions are defined by the new coupling constants associated to each of them. One curious outcome 
is that the larger the coupling (to be called from now on "flavour condensate" couplings) associated to the action 
terms for each flavour, the smaller the quark mass becomes. That is, the top quark shows the smaller condensate 
coupling. The modified Feynman rules determined by the new vertices are defined. The gluon free propagator remains 
unaltered, but the fermion ones decompose in two terms showing the same massive pole. One of them resembles a 
scalar particle propagator and the other shows a Dirac propagator structure but including the same massive pole in 
addition to the usual zero mass one. We postpone the study of the renormalization (the behaviour of the running 
strong and flavour condensate couplings) and start here to explore the first implications of the analysis by evaluating 
the gluon self-energy in the g 2 order in the gauge coupling and all order in the flavour condensate ones. The result 
for the gluon polarization operator became transversal as implied by the gauge invariance for this second order in 
the coupling calculation when all orders of the flavour couplings are included. Further, the transversal part of the 
polarization operator is employed to evaluate the vacuum potential (taking all the mean field values equal to zero) 
as a function of the flavour condensate couplings. The potential in this case, at variance with the one evaluated in 
[22j resulted to be bounded from below. Moreover, for small coupling values, in the here considered approximation, it 
shows a dynamical symmetry breaking generating non vanishing values for the quark masses. Since the renormalization 
group properties of the theory are delayed to a next study, it is not possible to make here clear cut phenomenological 
predictions. However, as a first exploration, the results for the potential were employed to select the g coupling 
and the dimensional regularization scale parameter /i in order to fix the mass for the top quark nearly to 173 GeV. 
That is, the flavour condensate couplings are suggested as playing a role similar to a Higgs field in the approach. In 
this calculation the values of the coupling g 2 and the scale parameter fj, result to be far from satisfying the known 
connection between the QCD running coupling and /i. However, it can be taken into account that in the considered 
two loop calculation, all the six quark contributions are "uncoupled". That is, all propagator entering a diagram 
have the same flavour. Therefore, the hierarchy of quark masses, if really implied by the theory, as we currently 
suspect, is to be expected only after including higher loop terms in the potential. In this case mixing among different 
quark species has the chance of determining a realistic dynamical symmetry breaking effect. This view about the 
need of including higher loop corrections coincides with the one exposed in Ref. [2l[ , where it was remarked that this 
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requirement becomes a natural one in a situation in which bound state effects are needed in determining a breaking 
of the flavour symmetry. In our considered case, to predict the quark mass hierarchy, it should occurs that a minimal 
effective potential appears for a large value of a single flavour quark condensate, while the other ones become very 
much smaller in size, and this minimum lays below another possible one in which various quarks condensate couplings 
get equal values. This situation only could have the chance of happening for the considered expansion, when the 
potential includes terms in which different sorts of quark propagators are present. Therefore, it is possible to conclude 
that the low values of the coupling got in the evaluation of the potential done here, is not yet reflecting a drawback 
for the analysis to be valid in describing the mass hierarchy. According to the above remarks, the next important 
steps in the search seems to be the study of the renormalization group properties of the expansion and the evaluation 
of three loops corrections to the potential for determining whether or not a quark mass hierarchy follows. 

The presentation proceeds as follows. In Section II the proposal of the local and renormalizable gauge theory 
including gluon and quark condensates and possibly being equivalent to massless QCD is presented. The case of 
retaining only the quark condensate is there shown to reduce to a modified local and renormalizable massless QCD 
action. Section III is devoted to evaluate the gluon polarization operator in second order g 2 of the strong coupling 
and all orders in the flavour coupling parameters defining the new vertices, and to check its transversality. Finally, 
Section IV making use of the transversal part of the polarization operator, presents the evaluation of the effective 
potential in the two loop approximation and all orders in the new interaction parameters. In the Summary the results 
are reviewed and possible extensions of the work commented. 



II. A LOCAL AND RENORMALIZABLE MODIFIED QCD 



Let us start by reviewing the main elements of the previous work which will be employed in the following discussion. 
The generating functional of massless QCD written in Euclidean variables as in Ref. [2(| has the expression 

7M nri ? 71- IlhJhMA 
z\j,v,v,U\- /[0;W ] ' 



dk . , . , 1 
(2nj 



ex P(/ 77^D](- k )oD{k)j{k))x 



^-p^ / (~fc)G / (fc)7 7/ (fc))x 



dk 

/ = 1,2,...6. (2) 



ex P( / 7^D^(~k)G gh (k)ak)), 



The functional is associated to the action S g depending on the gauge interaction coupling g. Sq for g — defines the 
free action. S g and the vertex part Lagrangian V mt are defined in terms of the six quark ^ j (/ = 1, ...,6), gluon A 
and ghost x fields in the form 

S g = [ dx{--F^ v - ±d,A1d v At+Tf^,Dfc b -J2 W f ^ D P% ( 3 ) 
J f 
V mt = S g -S , (4) 

where the field intensity, and covariant derivatives follow the conventions 

K» = ~ d »K ~ 9f abc AlAl : 

D« = + ig A^ , Df - d^ ab + gf abc A* , (5) 
{7„, lA = -28^, [T a T b ] = if abc T c . 



4 



The gluon and quark condensate parameters representing the modified vacuum free state used to construct the Wick 
expansion enter the generating functional only though the free quark and gluon propagators, which have the forms 

DfM = - (1 - a)^)9 N (k) + C g S D (k)S^, 

G«(jfe) = 0N ^ k \ +C f S D (k)I), 

G gh( k ) = 5 ~^2~- ( 6 ) 

The derivation considered the possibility of quarks already having the mass m appearing. These expressions are 
solutions of the tree approximation for the propagators, since the 4-dimensional Delta functions is a solution of the 
homogeneous Dirac and Klein-Gordon equations, because at zero momentum both equations vanish. As discussed in 
the Ref. , the functions 0/v(&) regularize all the propagators at zero value of the momentum as required by the 
gauge theory quantization [26| . The notation of Ref. [27| is closely followed along the work. The Einstein summation 
notation in the flavour indices is employed and a flavour (or other kind of index) between parenthesis means that it 
does not follow the Einstein convention. 

In Refs. [IB413; the alternative form of the functional integral (QJ for massless QCD was obtained after modifying 
the free vacuum, by creating in it gluons, quark and ghosts of nearly zero momentum. These changes of the initial 
condition in the process of connecting the interaction, after using the rules for constructing the Wick expansion 



( 28]), led to modified free propagators in ([5]). Some quantities had been evaluated with these modified expansion 

( 3 IS H3) and its general gauge invariance was argued 0, Q11 |2l|. It was possible to check that if the gluon 
condensate parameter defined here is chosen to match the estimated value of the more usual gluon condensation 
parameter < g 2 G 2 >, then the masses of the light quarks are predicted to be the constituent quark mass of nearly 
333 MeV. In addition sample loop calculations of the vacuum energy as a function of the condensate parameters 
were done in Refs. [18|, [lj, |2l| , in one loop order or by employing the ladder approximations for the quark and gluon 
propagators. They gave indications of a dynamical generation of the quark and gluon condensates. For gluons the 
one loop result gave a form of the potential similar to the potential in the early chromomagnetic Savvidi models of 
confinement. 

However, the working expression ([I]) of the generating functional shows two limitations, one of theoretical nature and 
another technical one. The technical one is associated with the fact that in the new contributions existed singularities 
requiring of a special regularization process to be well defined [l8|]. The theoretical one can be described as follows. 
In Ref. [20j we linearized the quadratic terms in the gluon and quark sources by means of introducing Gaussian 
integration over auxiliary fields. This led to a functional integral in which interesting dimensionally transmuted gluon 
and quark propagators already appeared. Next, in Ref. [13] we restricted the discussion to the case in which only the 
quark condensate parameter were retained, seeking to transform the modified Feynman expansion in a more helpful 
form by integrating over the auxiliary parameters. It was possible to transform the generating functional by mapping 
the effect of the condensate, from incorporating a Delta function term in the free propagator, to the addition of a new 
term in the massless QCD action associated to a vertex, showing two gluon and two quark legs. However, the form 
of this vertex, although being Lorentz and global color invariant was not local in space-time, since it had the form of 
a double integral over the product of two pairs of fields evaluated at two different space-time points. This effect was 
a direct consequence of the particular Dirac Delta structure in the momentum space of the free propagator. 

In this section we attempt to modified the construction of the modified expansion in order to overcome both of the 
problems. The starting idea for this purpose was suggested by the form of the generating functional obtained in Refs. 
|20l I22I ] . In beginning, let us consider the result for Z given in Ref. [2(| incorporating gluon as well as quark zero 
momentum condensates in the free ground state before connecting the interaction. The expression for Z is given by 
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Eq. (12) of Ref. [H in the form 

- Iff a a a a 
Vf,Vfid] = j? / dadxdx exp [ - ^ X/,«X/, U ^] 



exp| 



( i (-A ; )lD F (fc)i(A ; )+^r7 / (-fc)G F (A : )r ?/ (fc)]+^(-fc)G fl F h (A ; )C(A;)])], 



(2tt)^ 



X = (xi,X2, -Xf-iXe), X = (Xi,X2)-X/-)Xe)) ( 8 ) 
da^dx = II r d X>, r - (9) 

where the vertex interaction Lagrangian V OT * is defined in Q and in which the propagators for gluons, quarks and 
ghost fields are the usual Feynman ones as indicated by the superscript F. The whole expression is a Gaussian mean 
value over the gluon and quark auxiliary parameters a, x >X- They are constants in coordinate space, a has Lorentz 
and color gluon indices and x and X have spinor and quark color indices. Since the Delta function terms were linked 
with the zero momentum Fourier components of the sources, the number of auxiliary Gaussian variables required to 
represent in a linear form these zero momentum components of the sources is reduced to the small number of the 
indices of the internal field components. 

Before continuing, let us rewrite the expression for the vertex Lagrangian y mt appearing in the above relation 
in order to discuss the coming points. To simplify the writing let us define a notation absorbing the condensate 
coefficients in new defined parameters as follows 



X* = / 2C 9 ^ a gi , Cf a i f / Cf a-i 

p (2tt) d m ' 1,r (2tt) d f ' r ' f ' r (2tt) d 

Now, let us call S* , Sq the actions S g , Sq after being evaluated in the shifted fields which define Vmt in the formula 
(O through Vi n t — S* — Sq. These quantities, as expressed in terms of the fields and the above redefined condensate 
parameters, have the expressions 

^* ^ r a , 2C„ i .Cf.i — .Cfi 5 5 n 

= J dx[-\FZ v {A + 4>)F%{A + <t>) 

E (*/ n^( A + il^(AW f +f q i llM D^{A)%) 

f 

E $ *r„ ig4> a X i i f+ 
f 

W f z 7m i<^T^ + ^ i 7 „ i ff 0^ J '*/)]> ( 10 ) 

5 » = ^[^((^) i ^+( I ^) i x/,* + (^)^x / ,| ! ^] 



= I dx [ - \{d^A a v - d„A«)(d„A a u - 3 U AD - l-d^A^At - c a d 2 c< 



/ 

One important point in the above expressions is related with the last two terms in S*. Note that these are the 
only linear in the radiation fields terms appearing and that they vanish when the coupling is set to zero. Since the 
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condensate auxiliary parameters are constants in space-time, it follows that these terms are just proportional to the 
zero momentum Fourier components of the entering kind of radiation field. Therefore, when this terms are expressed in 
terms of the functional derivatives of the corresponding field sources, they are proportional to the functional derivative 
over the zero momentum component of the sources. Therefore, their action over the free generating functionals in 
([7]) vanish because the functionals are not depending of these components, thanks to the Nakanishi zero momentum 
regularization appearing factors in (flj . Since the condensate auxiliary parameters are constants they do not appear 
in Sq which due to the masslessness of the theory is defined by pure differential operators which action on them 
vanish. 

Henceforth, the action to be employed in what follows is 



S* g = S;[A,-f,^,c,c,a, X ,x} 

= I dx[-\F^(A + ct>)F- v (A- 



—d»Ald v Al ~ c a d»D a J>(A + 0)c b ) 



1 

2a' 



+ f f i^ig^K j l3j)}. (12) 

These last properties are useful to obtain an alternative functional integral form for Z. For this purpose consider 
the rewriting of this quantity as follows 

zb,wM= f[0>0|0>0] , 

_ Iff > a a a a 

V, % & C] = jf J J da dx d X exp [ - ^ x),uX).u ^1 x 

exp [-£o[— — - - — ]x 

5 j Sr] 1 —Srj — (5£ 

cxp[ / ^w~ fc )^( fc )j( fc )+E%(- fc ^ ( 13 ) 

Then, by returning each of the free generating functions to their functional integral original expression giving rise to 
them, it follows the relation 

eXp[ - 5 °^4'^4'^ ]X / m *'*'"' Cl 

exp[ / dx(C [A,^,^,c,^+j(x)A(x)+Y(rj f (x)^ f (x)+ : 9 f (x)r] f (x))). 
J f 

= J PL4, *,*,c,c] exp [ J dx(£o[A,*,%c,q-£o[A,*,*,c,q+j(x)A(x)+^ 

= J V[A,%%S,c]exp[ J dxY,(Vf(x)y f (x) + y f (x)ri f (x))}- (14) 
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Employing the above expression it follows for Z 



a a a a 



sri ^ s 5 5 5 5 8 _-, 

exp[S [—,—,— r ,—,—r-,a,x,X\\x (15) 
J oj or) —or) or) —or) 



J V[A,V,y,c,c]exp[ J dx (j{x)A(x)+J2(Vf(^f(x) + ^f(x)r) f (x))) 



= Jf J J dad xdx exp [ - X/, rX/, r " -y^] x 

/" P[A,*,*,c,c]exp[ J dx (S;[A,*,*,c,c,a,x,x}- 

]{x)A{x) +J2(Vf(x)^f(x) + *f(x)7 lf (x)))]. 
f 

where the action S* takes the form 

S* g = S*[A,%,y,c,c,a,x,x] 

r i on , on , 



i on ! 

f 



+ ^ 'K^xJ)]. (16) 

This formula is exactly the same that was obtained in Ref. [22l ] when the gluon condensate parameter vanishes. This 
can be seen more clearly after considering that all the derivatives of the quark auxiliary parameters vanish. This 
alternative expression suggests that the above defined local and renormalizable modification of QCD has the chance 
of being equivalent to massless QCD. For seeing this possibility it can be observed that the non locality and lack 
of gauge invariance of the theory defined by Z, are just determined by the coordinate independence of the auxiliary 
variables. In other words if x l f> X/ an< l a ^ are formally considered as arbitrary functions of the coordinates, the 
action S* (after substracting from it the gauge fixing and the ghost terms) becomes gauge invariant under the gauge 
transformations 

w{x)=exp[-i X a (x)T a }, 
w ad (x) = expH K(x)T% d ], (TZ d y k = i f iak , 

A'^x) = w ad A^(x)w^ - ^-w ad {x)d ll w ad (x)~ l , 

i&f(x) = w(x)^! f (x)w(xy 1 ,^ f (x) = w(x)~$ ' f(x)w(x)~ 1 , 

c'(x) = W ad (x)c(x)Wad(xy 1 ,c'(x) = W ad {x)c(x)w ad {xy 1 

<^(x) = w ad {x)(f)^(x)w ad (xy 1 , 
Xf = w{x)xf(x)w(x)~ 1 , 

Xf = w(x)x)(x)w(x)~ 1 , (17) 

where the auxiliary functions (before being the parameters) are represented as elements of the SU(3) algebra. The 
boson fields like A fi (x) — A^(x)T^ d , are expressed in terms of the adjoint representation for the algebra generators T® d , 
(the subscript ad indicates "adjoint") and the fermion quantities in the fundamental one in terms of the Gell-Mann 
matrices. 

However, the constant in the coordinates character of the condensate parameters, followed as a direct consequence of 
the particular nature in which the modified free vacuum state including condensed zero momentum gluons and quarks 
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was constructed. Thus, it is possible to conceive a wide range of alternative ways of defining generalized modified 
free vacuum wavefunctions within the manifold of states annihilated by the BRST charge of the free problem. Such a 
possibility, leads to the expectation about that a procedure can be developed in which the ending outcome can be of 
the similar form as (|15l) , but in which the Gaussian integrations over the condensate parameters become substituted 
by Gaussian functional integrations over condensate functions depending of the coordinates. We had attempted to 
determine such an alternative procedure, which validity will shows that the modified theory is in fact equivalent to 
massless QCD, but up to now we had not been able to find it. In this work however, we simply will propose the 
generalized expansion and consider the first implications of the local gauge invariant alternative forms of QCD that 
they represent. 

Therefore, we propose the following expression for the generating functional of a local and renormalizable gauge 
theory which, as it was argued above, has the chance of reflecting nonperturbative properties of massless QCD 

-Iff ^ a.f l (x)af l (x) 

Z\j,V,Tf,Z& = jf J J P[a,x,x]exp[-^x>, rW f< r (s) ~ M 3 ^ > 

J PL4,¥,*,c,c]exp[ J dx (S;[A,y,y,c,c,a,x,x} + 

j(x)A(x)+Y / (rJ f (x)^ f (x)+^ f (x) Vf (x)))], (18) 
/ 

in which the action is defined by 



S* = S*[A,i&,i&,c,c,a,x,X. 

1 1C i 1C 

dx[ -F^(A + {^B-FcftF^iA + (^) . 



- ^ d„A^d v A a v -cfd^iA + i— ^)>o») c") 



C 1 IC 1 i 

E ( frfnWt WV,]. (19) 



In the case of the vanishing quark condensate parameters, we expect that the expansion can lead to interesting physical 
consequences in the low energy region near 1 GeV, were a similar discussion in the previous non local expansion gave 
predictions for the constituent quark masses of the light fermions and reproduced the Savvidy chromomagnetic 
effective potential form, as a function of the gluon condensate parameter. Such results ca be expected to also come 
out from the generalized formulation proposed in the mean field approximation. 

However, in what rest of the paper we will explore the case being closer to the possibility for quark mass generation, 
by simplifying the discussion fixing to zero the gluon parameter. 

III. THE PURE QUARK CONDENSATE CASE 

When the gluon condensate parameter is set to zero, the Z functional takes the simpler form 

^(x)a a ^(x) 

/ 



J V[A,^,^,c,c]exp[ J dx (S* g [A,^,^,c,c, X ,x}- 



j(x)A(x) +Y / (Vf(x)^f(x) + ^ f (x)r, f (x)))}, (20) 
/ 
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where now the action gets the expression 

S; = S* g [A,*,*,c,c, X ,x] 

= Jdx [~\f; v F; v - ±- d ,Ald„Al-T*d,Dfc b ) 



- E (*} H^X j } ijgfc)* + (j§y>M (21) 

in which x'f an d Xf are now space-time coordinate dependent functions. Next, similarly as it was done in Ref. (2^ . 
we can make use of the fact that the dependence of the action S* is quadratic in the auxiliary functions Xf an d 
Xj, and thus the Gaussian integration can be explicitly evaluated by solving the Euler equations 

SS* [A, W, vt, c, c, y, y] C„ i „• 

St f = ~ X> " ( (2# )3 ^^*/ = °' (22) 

<55* c, c, y, y1 — 4-.. C„ i 

g - X/ + *'/ «7>;( (^U) 5 = 0, (23) 

B*=8i i d + igA%Ti i , = -6**5 + igAp* , (24) 

for the auxiliary functions Xf an d X/ an d substituting in (|20p . After doing this, the generating functional can be 
written as 

Z = ^ J V[A, % c, c, ] exp[5[A, *, ¥, c, c]] , (25) 
S[A,W,*,c,c] = 5 rogcd [A,W,*,c,c]+ (26) 
5 mgcd [A, c, c] = / ^(-If;^;, - - i 7 M^*/ ~ ^d fl Dfc b ), (27) 

where the action S mqc d is the usual one for massless QCD and as in Ref. (22[, new action terms appear, one for each 
quark flavour /. They have the expressions 

S*[A, = / d^j i 7p ^ J7 ,^ fc */' ■ (28) 

Note that we have slightly generalized the writing of the new action by a noting the gauge invariance and locality 
is not lost if we consider that the two quark fields appearing in the terms can have different flavours. However, 
if the matrix Cf /< is assumed to be hermitian, it can always be diagonalized by a unitary transformation in flavour 
space. Thus, we will assume that the matrix is diagonal. Therefore, it follows that in the case in which only the 
quark parameters are present, the resulting action is a modification of the massless QCD, in which new local and 
gauge invariant interaction terms appear, one for each quark flavour. The resulting theory has few properties that 
are underlined below: 

1) The new terms produce four legs vertices which curiously do not affect the renormalizability of the perturbativc 
expansion, as it might be expected due to the dimensional character of the new condensate coupling t^jit parameter. 
This property can be understood after noticing that the quadratic in the fields parts of the terms contribute to the 
new free quark propagator and made it to decrease with the square of the momenta in the ultraviolet region. Thus 
the power counting renormalizability of massless QCD is retained. 

2) In addition, these action terms seems to be possible counterterms appearing in the context of the renormaliza- 
tion of the standard massless QCD, possibly representing dimensional transmutation effects. Thus, the interesting 
possibility arises that a special renormalization procedure implementing a dimensional transmutation can also validate 
the theory proposed, as an effective action for massless QCD. The investigation of this question will be considered 
elsewhere. 

In what follows we will consider the evaluation of the one loop gluon self-energy and the vacuum energy as a 
function of the condensate parameters. 



IV. THE FEYNMAN EXPANSION 



Let us describe in this section the Feynman expansion for the case in which only the quark condensate parameters 
are non vanishing. For bookkeeping purposes we will consider this expansion in Minkowski space. Thus, let us define 
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an action in terms of the fields defined in Minkowski space, incorporating all the terms appearing in the expressions 
(|26l27l28p for the Euclidean action obtained in the previous section. Then, the Minkowski action to be considered in 
what follows is chosen in the form as 

S = [ I dx{-^F^F a ^ - ^-d^dvA™ + c a d ll D ab ^c h +Y^% *V-Dj?*J) 
- E j§jo I dx % ifo* l»D ik ^% (29) 
which follows the same conventions of Ref. [27]. The field intensity, and covariant derivatives are now defined by 

*%u = MS - d » A l - 9 f abc KAt, 

D» = d^ -igAp^, Df^d^ h ~ 9 r hc A% (30) 
{7m»7v} = 2s P v j [T a T b ]=i f ahc T c . 

The Minkowski action has been written with the same conventions employed in Ref. (27j and thus all the properties 
of the vertices and free propagators associated to the massless QCD part of the Lagrangian can be found in this 
reference. Then, let us only consider the definitions of the new form of the quark propagators and of the additional 
four legs vertices. Let us also take Cj ±, as Hermitian. Then, as mentioned before, by a global unitary transformation 
between the six quark fields, it can be brought to a diagonal form C/ f> = C/ St f . 
The expression for the quadratic in the quark field of flavour / Lagrangian is 



(2nY 

</..-"*, A, (d)¥ fl (31) 
Thus, the propagator Sf for each of the six quarks is given by the inverse of the Af operator appearing above. That 



is 



" - j^rnd 2 ) S f (x - y) = S(x - y), (>!2 . 

or in terms of their Fourier transform 



C 



{rP» + j0y5P 2 )Sj{ P )=I. (33) 



Therefore the quark propagator for each quark flavour / is given as 



V 'vr (2-k)°P > 



Sf(p) — c — c 

^_^-5<')( P ) + ^'«^ 



rrif m f 



(mi — p) (mj—p 2 ) p 

Of 



(2tt) 

It can be observed that the free quark propagators show poles of mass mf = c which are inversely proportional 



to the corresponding condensate parameter Cf. In addition, each fermion propagator decomposes in the sum of a 



"scalar" like term Si s \p), being equal to a massive scalar field propagator times the spinor identity matrix, and a 



( f) 

" fermion" like component Sf (p) having the same spinor structure as the Dirac propagator but also showing an 
additional pole at the same mass of the scalar like part. The expression for free generating functional for each quark 
flavour / takes the form 

Zo[V,V,f] =exp[ij j0rpr)f(p)Sf(p)Vf(p))- (35) 
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Let us consider now the vertices associated to the new higher than quadratic terms in the Lagrangian, after assuming 
flavour matrix C as not yet diagonalized 

C o f = E T^yiT / dx % 7T (7p^7^ + 7vfl"7^)*A + 
/1/2 

+ E 7§tf / dx % T a l l»A»X k l^i- (36) 
/1/2 

After transforming the integrals to momentum space the expression for the three and four legs vertices can be 
written as follows 

<?£2)((r 2 ,i 2 ,/ 2) ( fc i> k 2 , k 3 ) = g ^ipT^(-(k la rr s (^y^ + (rr s (k 2a rr 2 ), m 

1/(4) (v,a)(uA) , , 7 \ _ 2 W1/2 rpixirpii^/ ny lS ( v\sr2 (oo\ 

1/ (r 1 ,i 1 ,/ 1 )((r 2 , l2 ,/ 2 )^l' /C 2,K3,/«4j - 9 1 a 1 b Vl > Vf ) • \ M ) 

The diagrammatic representation of these vertices is illustrates in figure [TJ As usual, the momenta directions 
coinciding with the direction of the line associated to a quark propagator means evaluating it in the value of the 
momentum. When the direction is opposite the evaluation is in the negative of the same momentum. As remarked 
before the rest of the diagrammatic rules are the same ones as in Ref. [27j . The quark line, which is the only one 
changing its analytic expression in the discussion done here is also depicted. 




FIG. 1: The figure illustrates the diagram associated to the new quark gluon interaction vertices. The analytic expression of 
the three legs vertex is given in (1371) and the four legs one is defined by (|3S[) . 



V. ONE LOOP GLUON SELF-ENERGY EVALUATION 



As it was mentioned before in connection with the extension of the work, the most important next step seems to 
consider the renormalization properties of the proposed modified version of QCD. This issue needs a separate study to 
be done. However, let us consider here two evaluations that can be simply made finite in the Minimal Substraction 
scheme (MS), but without adopting yet definitive renormalization conditions in order to not disregard important 
elements not yet clarified. 

Specifically, we will evaluate the one gluon self energy in the second order in the color coupling and all orders in 

( 2tt) d 

the quark condensate parameters C/ = - m ' . The close related two loops contribution to the vacuum energy in 
the same approximation, will be also calculated. In what follows the presented results correspond with calculations 
using the Feynman rules of reference (27j for gluons and ghosts propagators and all the standard vertices of massless 
QCD, after complemented with the before defined rules for the new quark propagator and vertices. 

The diagrams associated to the mentioned approximation for the self energy and in which the quarks participate 
are illustrated in figure [5J The terms in which only gluons and ghost propagators and vertices enter are identical to 
the ones evaluated in Ref. (27j are not needed to the evaluated. In the one loop approximation under consideration 
the diagrams entering are sums of analytically identical terms for each of the flavours. Therefore, we will evaluate 
the contribution to the selfenergy given by only one flavour. The total contribution is the sum of the same expression 
obtained for one quark evaluated in the flavour couplings of each of the quark types. The expressions for all the 
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FIG. 2: The picture shows the five diagrams associated to the fermion contribution to the two loop gluon self-energy. They 
are evaluated in the text in the same order as they appear from left to right in the figure. 



appearing vertices and propagators associated to gluons and ghosts can be found in the exactly the same conventions 
used here in Ref. 27]. As an example, the wavy lines design the gluon propagator as given by. 



n ab (h \ _ 5 ab d tiu (k) 

d^u{k) = g^ v - (1 - a)^0^, 



(39) 



and the vertices showing filled dots are the standard quark gluon, ghost gluon and the two types of gluon vertices as 
defined in Ref. [27| . The two kinds of vertices determined by the new terms in the action are represented by two 
concentric open circles. In order to evaluate the polarization operator it is convenient to decompose this quantity in 
its transversal and longitudinal parts as follows 



IT^(fc) = 



D - 1 



k 2 

k^k v 



k 2 ' 



nf(fc) 



k^k i 



n£(A0. 



(40) 
(41) 
(42) 



It can be recalled that the fermion propagator (I34[) is the sum of one term being proportional to the spinor identity 
matrix and another one which is linear in the Dirac matrices. Also the new vertices both have a quadratic in the 
Dirac matrices structure. Therefore, it follows that the expression associated to the diagram G5 in figure ((2]), showing 
two of the new three legs vertices, decomposes in the sum of two expressions: one in which both propagators are of 
the ''''scalar" type and another in which both are of the "fermion" kind. The same happens for the diagram G4 in 
which two usual three legs vertices participate. The terms associated to G\ in which one of each of the new three 
legs vertex appears results in two times the expression in which one propagator is " scalar" like and the other one is 
of the "fermion" type. 



A. Diagram G5 

1. Scalar like propagators contribution ui^ b 

The analytic expression for graph in which both quark propagators are of the "scalar" type, after evaluating 

the color and spinor traces can be written in the form 



J D 2 

"6p„W- / (2 7 r)^(m2-(p-|)2)( m 2_( p+ |)2 ) - 
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which determines the transversal and longitudinal parts as 

5T { ' (d-i)J {2^i{m)-{p-iy){m)-{ P + iyy l44j 

U (s)ab (k) - -2a 2 S ab - ( A{j> - q)2 (45) 

2. Fermion like propagators contribution IlifJ^ b 

Writing the analytic expression and calculating the traces of the diagram G5 in which the quark lines are " fermion" 
like results in the expression 



1 2 X ab 2 f d P° (V + g 2 )(4p^ - q^qu + q 2 g^) - S{p.q) 2 g^ 



U U)ab (k) - ^--n 2 l) ah m 2 ( dp 
ll 5^ W - 2 9 f J (2tt) d 



(46) 



H (p _ £)2 (p + |)2 (m 2 _ {p |)2)(m 2 _ {p+ ^ ■ 

In this case the integrals for the transversal and longitudinal functions become 

U)ab 1 g 2 5« b f dp D (4p 2 + q 2 Mp 2 -^) + g 2 (D-l))-8(p. q ) 2 (D-l) 

5T W 2(25-1)7(2^ (p - l)2 {p + 1)2^2 _ {p _ q _ ) 2 ){m 2_ {p+q _ ) 2 ) > i 47 ) 

n(/)« b (fc) = / 4(V-^x^ 

5i lJ 2^° g 2 _/ ( 2 7r)^(p-|) 2 (p+|) 2 (m2 -(p-|)2)( m 2_( p+ |)2)- ^ 

The above integrals which determine the contribution to the diagram G 5 to the gluon selfenergy, in general contain 
more momentum dependent factors in the denominators that the ones in massless QCD at the same one loop level. 
However, those integrals, and the ones appearing in the rest of the contributions to the self-energy, can be systemati- 
cally reduced to linear combinations of one loop scalar integrals, by employing the following definitions for the factors 
in the denominator determining the poles of the integrands 

D, = {p-\f, D 2 = ( P + q -) 2 , 

D 3 = (m 2 - (p - §) 2 ), D A = (m 2 (p + §) 2 ). (49) 



These relations can be inverted to express the quantities p 2 and p.q in the numerators as linear functions of q 2 and 
the D terms in various ways as follows 

P 2 = m 2 -^- 1(£> 3 + £> 4 ), p.q = \{D 3 - 1)4), 

m 2 f n 2 1 m 2 f 1 

p 2 = ^-- q ---{D i -D 1 ), p.q=-L--( Dl+ D i ), 

o m 2 f a 2 1 m 2 f 1 

P =-£-~\- 2 {Ds ~ ° 2) > p - q = ~^t + 2 {D2 + ° 3) > 

p 2 = -^+ l -{D l +D 2 ) 1 p.q= l -(D 2 -D 1 ). (50) 

Then, the following identity makes the work of decomposing the above Feynman integrals in simpler ones 



(p - § ) 2 {P + I) 2 (m 2 - (p - f ) 2 )(m 2 - (p + f ) 2 ) 

= m ) 
D1D2D3D4 

1 1 1 1 

~ D X D 2 + D 3 D 4 + D1D4 + D 2 D 3 ' ( ' 
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Since all the integrands in the transversal and longitudinal parts are functions of p 2 and p.q, all of them can be 
expressed as functions of the D factors and the square of the external momentum q, The following relations between 
the basic integrals resulting in the various evaluations done here follow 







dp D 


Di 


(2ir) D i 


D 3 


dp D 


D 2 


{2%) D i 


D 3 


dp D 


D 3 


(2it) d i 


Di 


dp D 


D 2 
u 3 


{2n) D i 


D A 


dp D 


D\ 


{2rr) D i 


D A 



dp D 


1 


(2n) D i 


D 2 


dp D 


1 


(2n) D i 


D 3 


7 D 

dp u 


D 2 


(2it) u i 


Da 


j D 

dp 


Di 


{2n) u i 


Da 


dp D 


Da 


(2ir) D i 


D 3 


dp D 


D\ 


(2ir) D i 


D 3 


dp D 


D\ 


(2ir) D i 


D 3 



dp 1 



(2tt) d iD 2 -' 



L ^)=l7^r^ = j^l^ (52) 

m 2 L m (m f ), 
(m 2 f +q 2 )L m (m f ), 
-q 2 L m (m f ), 



.2,2 , ., "'/'/ 2 



((q 2 Y + 4^-)L m (m f ), 



HI 



2 2 
f q 



((mj + q y + A^-)L m (m f ). 

Employing the results in references [27|, [29|, |3Q|]j performing the Wick rotation and analytically integrating the 
appearing Feynman parametric integrals with the use of the Wolfram Mathematica code, the three basic integrals 
appearing in the above formulae can be evaluated. The tadpole integral has the form [23, [29j 

f dp D 1 mf~ 2 B 

= J ^gTFf) = (i^ r(1 - (53) 

The massive scalar self-energy after making use of the formulae in Ref. [27], performing the Wick rotation and 
again analytically integrating the appearing Feynman parametric integral can be written as follows 



L 34 {q, m f ) 



dp 1 



(2n)° {m 2 +{p-l) 2 ){{m 2 +{p+l) 2 ) 

-Vr(2 - ^-)mf~ 2 [ dx(l + x(l x)^- 2 
(4tt)tt 2 J J mj 

T(e) , m 2 y 4m / + g\ , 
> ) 



(4ir) 2 - e Am 2 f + q 2 q 

(B[\(l+ . 1 ),l-e,l-e] 
/ Am r. + q 2 



2 //Irn2 



B[Ul- , q ),l-e,l-e]), (54) 



n'2 



+ 1 

where B[z, a, b] is the Incomplete Gamma Function and as usual e = 2 - 

B[z,a,b]= [ dtt a - 1 (l-t) b - 1 . (55) 



D 
2 ' 



Next, the self-energy term including one massive and one massless scalar, was also evaluated by employing the 
formula given in Ref. 30] by after the Wick rotation, analytically integrating the appearing parametric integral. 
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The result becomes 

dp 



Lu(q,m,f) 



(2 7 r)^(p-|)2(m 2 / -(p+|) 2 ) 

2 



-> ' r(e)r(i-e). , , „ 2 ,_ 2e g 2 



(47r) 4 1 (2 — e) J to^ — 

in which 2-Fi(a, c, 2) is the Hypergeometric Function 



(56) 



z? < h \ ^ (a)k{b)k z k T(a + k) 

2 F 1 (a,6,c,^)=^-^— -, (a) fc = ^-. (57) 



The massless scalar self-energy integral results in 



1,12 (?) / (2n) Di <•• 1 '-> " - 1 ' 2 



(27r)^(p-|) 2 (p+f) 2 

^-v-t^r ^-y . (58) 



Finally, the various contributions to the polarization operator associated to the diagram G5, for each type of quark 
flavour /, can be written as explicit functions of the quark masses rrif (or the quark condensate parameter C/) and 
the space time dimension D as follows. The amplitudes of the transversal components result in the form 

q 2 S ab m 2 

H^ ab (q) = ^ ' {D _ 1 J ) [(*rn 2 f + (D- 2)q 2 )L 3i (g, m f ) - 2L m (m f )] , (59) 
nit\<l)=^ 2 S ah m}L m (m f ), (60) 

n^ T )a6 (g) = T$(q) + T 3 i(q) + TfJ(g) + 1$(q), (61) 
T$(q) = 0, (62) 



a 2 S ab ( 2 \ 
n h M = - 2 {D-1) [ H D - 5) L m(m f )+^mj + (D - 2)q 2 )L 3i (q,m f )j , (63) 

T${q) = T£{q) = - J ,J (2m}(3 -D)-2-f + 2q 2 (D - 2))L u (q,m f )+ (64) 



a 2 TO f TO f a 2 

((l+-^r)(6 - 2D + 2-f ) + 4(£> + 1) + 2-f + 2{D -2)\ 
mj q z q z mj 

\{m 2 + 2(1 + hq 2 + ^)L m (m f )}, 
q z 1 D m~j: 



where the indices (1, 2), (3, 4), (1, 4) and (2, 3) here and below will indicate the contributions associated to the corre- 
sponding four terms in the decomposition (|5"Tj) of the factor DlD ^ DaDi ■ 
The longitudinal components get the expressions 

(65) 
(66) 

(67) 
(68) 



-£ ab (k) = 


U$(q) + 




U$(q), 


U$(q) = 


0, 






Ug(q) = 


g 2 S ab 4(l 


2 

+ ^) L m (m f ), 




U?hq) = 


U$(q) = 


-2g^ b [{^{m) 




(- 




+ l))L m (m f )]. 
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B. Diagram G2 

This is the simplest of the calculations. After writing the analytic expression of the graph by following the Feynman 
rules and evaluating the color and spinor traces, the polarization operator contribution becomes 

Tjab 1 \ a 2 <rab uu f dp D 1 

n^(«) = -4ff* f J {27:)Dim} _ p2 

= -4g 2 S ab L m (m f )g^, (69) 
n§r(9) - n5i(fc) - -4g 2 5 ab L m (m f ), (70) 

in which the transversal and longitudinal parts becomes equal and proportional to the scalar massive tadpole integral. 



C. Diagram Gi 

In this case, as noted before, since there is one new three legs vertex in the diagram, in which a product of two 
Dirac gamma matrices enter and another of the usual massless QCD which only contains one, the calculation reduces 
to two times the one in which one fcrmion like propagator and one scalar of the scalar type are employed in the 
internal lines. Following the same steps as before, the integral giving the self-energy contribution associated to the 
quark flavour / takes the form 

LL llu ,W-2gd m f J (2?r)D . (p + f)2(m 2 _ (p _ |)2)(m 2 _ (p + |)2) ■ (71) 

After getting from it the transversal part and expressing the integrand in terms of the D functions and q 2 , allow to 
obtain the analytic expression of this transversal part in terms of the above given scalar integrals in the form written 
below 

7 SZ(l) + W?Z(q) + W2£(q)), (72) 

(73) 

)-2) q 2 + 4m})L 34 {q,m f ) + (^ -6) L m (m f )j , (74) 
i f ), (75) 



Uf T (q) 


= m'iq) 


W$(q) 


= o, 


Wf 4 {q) 


g 2 S ab 
(D-l) 




8g 2 S ab 

~ D L 


W$(q) 


W£{q) 


2g 2 8 ab 
~ (D-l) 1 



— [(mj(3 -D)+ q 2 (D - 2) - ^l)L 14 (q,m f )+ (76) 

(Dl 2 + m)) T 

-2 L m (m f )\. 

Following the same steps as before, the evaluated expression for the coefficient of the longitudinal component 
becomes 

n?i(g) = - (X$(q) + XfM + Xf 4 (q) + X$(q)), (77) 
X$(q) = 0, (78) 

Xf 4 {q) = Ag 2 S ab (l + |) L m (m f ), (79) 

*i a l(<Z) = -^imK), (80) 

9 2 
TTW 777- 

X£(q) = -2g 2 5 ab [-f(m 2 q 2 )L 14 (q,m f ) + (2 - -^)L m {m f )]. (81) 

The contribution of the diagram G3 exactly coincides with the just evaluated for G\ . This can be seen by performing 
a mood change of the 

momentum integration variable in the expression for Gi. 
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D. Diagram G4 



The process of evaluation of the diagram G4 follows the same steps as the ones for G5. The form of the results 
for the integral after the traces are calculated is 

n 4 t,(fc) = n^ ab (fc) + nS a "(fc), (82) 

lC a6 (fc) = -2 9 2 S ab m}g^ I (2 f )D Z (m 2 _ (p _ | )2)(m 2 _ {p+ a )2) ( 83 ) 
= -2g 2 5 ah m 2 } L 3i {q,m f )g fJtu , 

n (B)a6,, x = _ /■ JP^ gp^ - 2q^ + (g 2 - 4p 2 )g„„ 

4 ^ 1 J 5 2 7 (2tt)^ (p - + |)2( m 2 _ (p _ |)2 )(m 2 _ (p + | )2) ' 

and the explicit formulae for the coefficients of the transversal and longitudinal parts of II^, 6 are 

Ii { § ab {k) = n { 4 i )ab (k) = -2g 2 S ab m 2 L 34 (q,m f ). (85) 
Similarly, the result for the transversal part of the contribution n^ ab in terms of the basic integrals become 
nft(«) = Bf 2 (q) + BfM + B$(q) + B? 3 (q), (86) 
Bt b (q) = - g ^g-^8 ab q 2 L 12 { qi m f ), (87) 



(D-l) 
9 



Q 2 S ab m 2 f „ „ 2m 4 f 

<(<z) = B 2 1(5) - - 2 p_ 1} [(2(2 - 73) + 2^)L m {m f ) + (m 2 (3 - 73) + 2g 2 (73 - 2) - -^)L 14 (q, m f )}. (89) 

Finally, the longitudinal term get the expression 

n&(«) = Cf 2 (q) + Cf 4 (q) + CfM + Cf 3 (q), (90) 

C$(q)=0, (91) 

C&(q)=2g 2 6 ab m} L 3i (q,m f ), (92) 
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CfM = Cf 3 (q) = -g 2 S ab m 2 (-f - l)L 14 (q,m f )+g 2 S ab -fL m (m f )}. (93) 



E. Gauge invariance and total transversal polarization operator 

After adding all the above evaluated contributions to the longitudinal component of the selfenergy it follows 

Hf («) = Uf L (q) + Uf L (q) + Uf L (q) + Ylf L (q) + Hf L (q) 

= 0. (94) 

Therefore, the transversality property is satisfied by the self-energy as it should be since it is a Ward identity which 
must be satisfied by any correction to the self-energy given by a well defined order in the perturbation theory. Since, 
here we evaluated all the terms in the expansion which are of order two in the strong coupling being and exact (to 
all orders ) in the flavour condensate parameters, the transversality should be satisfied. 

Further, after also adding all the transversal contributions and simplifying the result, the total unrenormalized 
one loop gluon self-energy can be written in the form 

U; b v (q)=ILf(q)(g^-^), (95) 

nf(q) = nf T ( q ) + nf T ( q ) + nf T (q) + n&(«) + n&( g ) (96) 

- J (73(73 -26) + 8)L m (m f ) + 9 -± - ) q 2 L 12 (q) + 



73(73- 1) K ' ' ~>~"^->> (73-1) 

+ 2{D-\) ^ W{B + ll)m2 f ~ 9D{2 + D ^ L ^ m /)' 
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Employing the Minimal Substraction renormalization scheme (MS), the one loop divergences can be eliminated 
by deleting the terms showing poles in the expansion of the above expression in a Laurent scries in the parameter 
e = ^p-- However, the resulting expression is cumbersome and since we will not analyze it here, is not written. After 
adding to IIy b (g) the free part of the gluon equation of motion, the following expressions for the gluon one loop 
propagator and its inverse follow 

Df-\ q ) = ^V(v - (i - a)3e£) - {q){g» v - (97) 



= tf-W r \q)){g„ v -^)+a^ 2 . 

H#(q) = 5 ab n T (q), (98) 

D ab (a) = — (a _ Ml) + —Ml fog) 

U » M} (q* _ n T (g)) q 2 )+ a q* ' [ ™> 

VI. VACUUM ENERGY 

Let us consider below the main steps in the evaluation of the one and two loops corrections to the effective action 
which will give the negative of the vacuum energy as a function of the quark masses to/ in that approximation. 

A. One loop term 

The one loop term reduces to the logarithm of the determinant of the new form of the inverse quark propagator, 
which is the result of the functional integral defining free quark generating functional. The calculation of this term, 
following usual steps in this special case, is sketched below 



T^{m f ) = i^logdetG- 1 ^) 



i 
p 



d l>" i.,„j...r '"I „ 



-^E^ogdet^a 



m f - p2 



TrD ST^ d P° , , r m f n VD V- d P° i , rP 2_m /n 

= -V E (2^7 ^ de '[^'] - — E (2^J logdet!^/] 

- v E <^ >os'K - ? 2 i + " D E (ioo) 



After taking the derivative the expression over m 2 , it follows 



= V D 2N L m {m f ) 

n r(i - #)to^- 2 

= ^ D 2iv 2 V , (ioi) 

(47r)~ 

which when integrated over m'j in the interval (0, m 2 ) by assuming e being in a region in which its real part is negative, 
leads to 

r (0) K) - V D 2N ) '\ f \ ■ (102) 

(47r)^(2 - e) 
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The finite part of the above expression in the Laurent series expansion can be evaluated by employing the expression 
of V D in terms of the four-dimensional volume and the dimensional regularization scale parameter jj, 



V D = ii A -°V A = (103) 



Then, the finite part the one loop effective action in the MS scheme becomes 

r £K) = 2 ^ lo ^) - (3 - 7 + log(4*)) 

= -Vj&imf). (104) 

The result after changed its sign, gives the dependence on mf of the vacuum energy in this approximation. The 
potential decreases at large mass values indicating the dynamical gene ration of the mass in this approximation. This 



outcome for the first correction was also obtained in Refs. |19l |20L l22| 



B. Two loop terms 

The two loop vacuum energy can be readily evaluated starting from the already known transversal part. This 
follows, thanks to the formula 



r (2) (m/) 



(D-l) f dq D Wj?(q) 



{2n) D i q 2 
(D-l)(N 2 -l) f dq D U T {q) 
{2n) D i~tf~' 



(105) 



which follows after noting that the diagrams defining the two loop approximation for the effective action can be 
organized as the sum of the diagrams defining the self-energy contracted with the gluon free propagator. 
After substituting in the above formula the evaluated expression for LTt it follows 

(2) g 2 (iV 2 -l) , mr> f dq D L m (m f ) g 2 5 ab (2 - D) f dq D 

T ( V/) = ^—(D(D 26) + 8) J g2 + ^-^ J j^MM+ 

n 2 (N 2 -l) f dq D 2D(D + 17)m 2 
+ 4 JWWi { g2 f +9D{2-D))L 34 (q,m f ). (106) 

In the above formula, the first and the second terms vanish in dimensional regularization. The third one can be 
expressed 
as follows 

rW( m/)= | 9^(N 2 - l)D(2 - D)V D f dq D dp D 1 



(2ir) D i (2n) D i (m 2 - q 2 )(m 2 - p 2 ) 



g 2 (N 2 -\)D{D + n)m 2 f dq D dp D 1 



{2n) D i (2n) D i q 2 (m 2 - q 2 ){m 2 - (p - q) 2 ) ' 



(107) 



Each of the two integrals appearing in the first term is equal to the L m in (1521) . The second term is proportional 
to a two loop scalar Master Integral given in Ref. [2{J. Then, the unrenormalized two loops effective action has the 



form 



r(2)(m/) = + 9(N*-l)D(2-D)V°ry-J) 
K " 4 (47r)-° f 

(N 2 -l)D(D + 17)(D-2)T 2 (2 + e)T 2 (l~§) 2 2D _, D 



2 2 2 W(D-3) y "V 

Substituting the formulae for the coupling g and the D dimensional volume according to 

g = go» D , V D = ^V\ 



g*m?-*V u . (108) 
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it follows 



2 2D-4irD 2 T /4 4 -4e 

g m f V = g V m f m f . 



Afterwards, expanding the formula (|108[) in Laurent series around e = 0, deleting the pole terms and tending to the 
limit e — >■ leads to the finite value of the two loop contribution to the effective action in the MS scheme 



r( 2 W) = -273.18 

V 31 47T 1 



322.47 



An 



A -r , m f \ 

m}Log(-jj-)- 



132.527 m 4 ^ Log 2 (^' 
; 4tt /i ' 

The total effective potential is then given by the expression 



V(m/) = 0.0656145 



273.18 m 4 f ^2- - (0.0379954 

/ 47T 



, m 



4 <?0 



322.47 2iL) m\ Log(^) + 132.527 m 4 ^ Log 2 (^) 



An 



4tt 



(109) 
(110) 

(111) 
(112) 



In Figure |3] the potential divided by the fourth power of the renormalization parameter \i is plotted as a function of 
the mass divided by /i for various small values of the coupling constant g . For higher values the coupling minimum of 
the minimum of the potential tends to disappear. However, as it was remarked before, the two loop approximation is 
insufficient to define whether or not the scheme predicts a hierarchical flavour dynamic symmetry breaking. Therefore, 
the evaluation done should not be expected to be relevant for answering the main physical question: the possibility of 
large quark mass generation at the values of the strong coupling. However, the results shown in Figure [3] indicate 
that at small values of the gauge coupling, the interaction is able to generate large masses for the fermions. The figure 
corresponds to fix as an example fi = 1 GeV , then it follows that for a small coupling near the value go = 0.0271828 
the fermion mass gets a value near the top quark mass m/ = 173 GeV. This evaluation suggests the interesting 
possibility that, if the theory is in fact equivalent to the massless one, the generation of large mass mechanism could 
also work for the low coupling electroweak scale. 




-6xl0 6 t 



FIG. 3: The figure illustrates the effective potential divided by ^i 4 as function of the ratio At /i = 1 GeV various small 
values of the coupling constant around go = 0.0271828 were chosen to evidence that the minimum of the potential can be fixed 
at a m/ being close to the top quark mass of 173 GeV. For higher values of the coupling the minimum tends to disappear. 
However, as described in the text, the lack of the possibility for generating a hierarchy of mass in the considered two loops 
approximation, determines the need for improving the evaluation of the potential, in order to decide whether or not a large 
and hierarchical mass generation effect can be predicted at the large values of the strong coupling observed in Nature. 



VII. SUMMARY 



We have proposed an improved version of the modified version of QCD discussed in previous works, which shows 
local gauge invariance and include the same kind of gluon and fermion condensate parameters. The analysis done 
for constructing the proposal suggests its equivalence with massless QCD, after considered with a special renormal- 
ization procedure designed to implement the dimensional transmutation effect. In respect to the gluon condensation 
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properties, the theory should reproduce the previous derivation of the constituent mass for the light quarks given 
in [I?], [H| and the Savvidy kind of potential as a function of the gluon condensate parameter. The study of the 
possible improvement in the predictions determined by the gained gauge invariance and locality of the description 
will be considered elsewhere. For the case of only retaining the fermion condensate parameter, the fermion auxiliary 
functions were integrated leading to a theory described by a new Lagrangian given by the massless QCD one but 
including a new gauge invariant term for each quark flavour. The new Lagrangian terms can be considered as local 
modifications of the ones derived in Ref. [22] . In the same form, they have a similar component constituted by the 
product of two gluon and two quark fields, but now evaluated in the same space-time point. The new terms also 
show components leading to three legs vertices which were absent in the previous form of the generating functional. 
The terms determines masses for all the six quarks which are given by the reciprocal of the new Savour condensate 
couplings linked with each quark type. Therefore, the strength of the condensate couplings decreases with the masses 
of the associated quarks. The gluon self-energy was evaluated up to the second order in the gauge coupling including 
all orders of the flavour condensate ones. The result, being associated to a well defined order in the couplings, satisfies 
the transversality condition as required by the gauge invariance. In addition, it is also gauge parameter independent. 
The transversal part of the self energy is employed to evaluate the two loop contribution to the effective action at 
zero mean fields (minus the vacuum energy) as a function of the flavour condensate couplings. The transversality and 
gauge parameter independence of the gluon self-energy, then also determined the gauge parameter independence of 
the evaluated potential. Afterwards, in order to interpret the result, we assumed that the introduced flavour couplings 
represent dynamically generated quantities, as it was effectively the case in the precedent works for the similar non 
gauge invariant couplings which motivated their consideration in this study. Then, in this first approximation, the 
potential is able to predict the dynamic generation of quark masses, for sufficiently small values of the QCD coupling. 
It can be remarked that the generation happening only for small couplings at the considered approximation, does not 
mean a negative answer to the question about whether or not the scheme can address the quark mass hierarchy. This 
is can be understood by noting that in the two loop order, the diagrams can not yet incorporate lines associated to 
two different kinds of flavours. For this appearance to happens, at least three loops corrections are required. It is 
reasonable that in order to explain the quark mass hierarchy as a dynamic flavour symmetry breaking, there should 
exist "interference" like effects in the vacuum energy corrections. In them, the contributions of diagrams showing 
two or more kinds of fermion lines might tend to rise the energy of the configurations with equal values of the quark 
condensates, making them more energetic that the ones in which one of the quark condensate parameters gets a large 
value and the others hierarchical lower ones. At the moment we have the impression that the considered framework 
is ideal to realize the Democratic Symmetry Breaking properties of the mass hierarchy remarked in Refs. 0,3- The 
predictions of the general discussion including gluon condensates for the low energy processes, as well as the renormal- 
ization properties and the evaluation of three loop contributions in the case of only having the fermion condensate, 
are expected to be considered elsewhere. 
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